Static Aeroelastic and Longitudinal Trim Model of Flexible 
Wing Aircraft Using Finite-Element Vortex-Lattice Coupled 

Solution 

Eric Ting* * 

Stinger Ghaffarian Technologies, Inc., Moffett Field, CA 94035 

Nhan Nguyen^ 

NASA Ames Research Center, Moffett Field, CA 94035 
Khanh TrinlE 

Stinger Ghaffarian Technologies, Inc., Moffett Field, CA 94035 


This paper presents a static aeroelastic model and longitudinal trim model for the analysis of a flexible 
wing transport aircraft. The static aeroelastic model is built using a structural model based on finite-element 
modeling and coupled to an aerodynamic model that uses vortex-lattice solution. An automatic geometry 
generation tool is used to close the loop between the structural and aerodynamic models. The aeroelastic 
model is extended for the development of a three degree-of-freedom longitudinal trim model for an aircraft 
with flexible wings. The resulting flexible aircraft longitudinal trim model is used to simultaneously compute 
the static aeroelastic shape for the aircraft model and the longitudinal state inputs to maintain an aircraft trim 
state. The framework is applied to an aircraft model based on the NASA Generic Transport Model (GTM) with 
wing structures allowed to flexibly deformed referred to as the Elastically Shaped Aircraft Concept (ESAC). 
The ESAC wing mass and stiffness properties are based on a baseline “stiff” values representative of current 
generation transport aircraft. 


I. Introduction 

In recent years, considerable attention within the aircraft design industry has been directed at incorporating lighter 
weight materials in the construction of aircraft structures. These efforts target reduction of aircraft weight, which 
translates into a lower lift requirement and in turn, reduces induced drag and thrust requirements. Usage of modern 
light-weight materials such as advanced composites has been adopted for use in the design of new aircraft. These 
materials are able to provide weight savings while maintaining the same load-carrying capacity as older material 
selections, also allowing the structural rigidity of the designs to be reduced. It becomes increasingly important for 
these modern designs to take into account the aeroelastic interactions between flight aerodynamics and the flexible 
aircraft structures within flight. Understanding and modeling these interactions can aid engineers in the analysis of 
design selections incorporating these lighter weight materials. 

In 2010, a conceptual study titled “Elastically Shaped Future Air Vehicle Concept” 1 was conducted to investigate 
the benefits of several advanced aircraft concepts over a conventional design. The study showed that there exists 
potential benefits in shaping wing surface aeroelastic deformation actively in flight with active control. The possibility 
of increasing aerodynamic efficiency by using active wing shaping for modern wing structures with reduced structural 
flexibility was realized. Development of modeling tools that can be used to investigate flexible wing aircraft and 
their aeroelastic behavior is emphasized in order to pursue these areas of research. These tools can also lead the way 
to flexible wing design optimization and the development of novel control surfaces to achieve active wing shaping, 
such as the Variable Camber Continuous Trailing Edge Flap system being investigated in a joint effort by NASA and 
Boeing. 2-3 
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This paper outlines the development of a static aeroelasticity model and a three degree-of-freedom longitudinal trim 
model for an aircraft with flexible wing structures. The static aeroelasticity model utilizes a one-dimensional structural 
model of the the wing structure as a beam in coupled bending -torsion. 4 5 The aeroelastic model also takes into account 
engine thrust forces and the effect of aero-propulsive-elasticity. 5 Previous studies have utilized the Galerkin method 4 - 6 
to formulate a discretized weak-form solution to the structural equations. The similar numerical technique of finite- 
element method (FEM) 7-9 utilizing shape functions will be used in this study, as approached by previous work. 5 An 
aeroelastic model is generated by coupling an aerodynamic model based off vortex-lattice data with the structural 
model through a geometry modeling tool that can create aeroelastically deformed aircraft models. 

A three degree-of-freedom longitudinal trim model is developed as an extension of the static aeroelastic model. 
The aircraft’s trim state is determined by balancing the longitudinal forces and moments on the flexible wing aircraft 
such that equilibrium is obtained. The trim model uses a standard Newton’s method approach to solve the aircraft’s 
equilibrium equations in the aircraft stability axes, while accounting for aeroelastic deflections. The modeling capabil- 
ity allows for the longitudinal states and the wing aeroelastic shape at trim to be determined. A comparison between a 
rigid wing aircraft and the effect of incorporating aeroelastic effects can be conducted. 

The static aeroelasticity and the three degree-of-freedom longitudinal trim model are applied to various different 
aircraft models. The models are based on the NASA Generic Transport Model (GTM), which is a transport aircraft 
airframe similar in class to the Boeing 757. 10 A simplified model based upon the planform and the characteristics of the 
jig-shape wing is developed, called the “Idealized Wing Alone” model. Validation of the developed static aeroelastic 
framework is conducted using the Idealized Wing Alone model and the results are compared against NASTRAN 
aeroelastic results provided by Boeing Research and Technology. 2 A full aircraft model with fuselage, tails, and 
engines referred to as the Elastically Shaped Aircraft Concept (ESAC) is utilized. The ESAC model can utilize a rigid 
wing model or wing model where aeroelastic deformation is considered. For the flexible wing model, the ESAC model 
utilizes mass and stiffness values representative of the GTM wing. The static aeroelastic framework is applied to the 
ESAC, and the resulting model is extended to develop the trim solution capable of handling wing aeroelasticity. 

II. Aircraft Models 


A. Elastically Shaped Aircraft Concept 

The Elastically Shaped Aircraft Concept (ESAC) is a model of a complete, typical, transport aircraft configuration. 
It is developed based upon the NASA Generic Transport Model (GTM), 10 which represents a notional single-aisle, 
mid-size, 200-passenger aircraft. The GTM is a research platform that includes a wind tunnel model and a remotely 
piloted vehicle. An extensive wind tunnel aerodynamic database 10 also exists for the GTM configuration lending 
itself for understanding and validation of results. Figure 1 is an illustration of the GTM planform. The benchmark 
configuration represents one of the most common types of transport aircraft in the commercial aviation sector that 
provides short-to-medium range, 3000 nautical mile, passenger carrying capacities. 



Figure 1. Benchmark GTM Planform 
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The geometry of the ES AC is obtained by scaling up the geometry of GTM wind tunnel model by a scale of 200: 1 1 . 
In the aeroelastic model, the jig-shape wings of the ES AC can be analyzed using a rigid model or be allowed to freely 
deform based on reference wing stiffness values. The equivalent beam model for the wing and the aircraft mass 
and baseline stiffness values are built using a component-based approach. The aircraft is divided into the following 
components: fuselage, wings, horizontal tails, vertical tail, engines, operational empty weight (or OEW equipment), 
and typical load including passengers, cargo, and fuel. The fuselage, wings, horizontal tails, and vertical tail are 
modeled as shell structures with constant wall thicknesses . 11 

B. Idealized Wing Alone Planform Model 

Development of the Idealized Wing Alone model is conducted using the undeformed jig-shape wing of the ESAC. 
Using the geometric wing pre-twist, measurements of mean camber line, and planform shape, the removal of the 
fuselage, tails, engines, and pylons allows for the construction of the wing alone model. This wing alone model is 
idealized as a surface with no thickness, and the wing root is extended from the planform of the original wing to the 
aircraft centerline. The idealized surface can then be twisted based on the jig-shape geometry wing pre-twist, or shaped 
based on camber line measurements. The resulting Idealized Wing Alone model represents only a wing surface, and 
the jig-shape’s inherent twist and camber can be individually activated on the model for incremental analysis of their 
effect on static aeroelasticity. 

Figure 2 is a depiction of the Idealized Wing Alone model within the vortex-lattice modeling program Vorview. 



Figure 2. Idealized Wing Along Planform In Vorview 


III. Wing Structural Modeling 

A structural model of the wing using beam theory is developed and is later incorporated into a fully coupled 
structural-aerodynamic aeroelasticity model. The model is similar to those developed previously in other studies . 4,5 

A. Reference Frames 

Figure 3 illustrates three orthogonal views for a typical transport aircraft and several associated reference frames. 
These reference frames are useful in developing the structural models of the lifting surfaces of an aircraft, although 
the coordinate frames associated with the aircraft wings are primarily used in this analysis. The aircraft body-fixed 
reference frame B is defined by the unit vectors b \ , b 2 , and £ 3 , which are aligned with the aircraft roll, pitch, and yaw 
axes, respectively. 
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Figure 3. Aircraft Reference Frames 


The reference frame C is aligned with the right wing’s elastic axis and is defined by the unit vectors c i, C2, and 
C 3 . Let A be the sweep of the wing’s elastic axis. The B frame can be related to C through three successive rotations: 
1) the first rotation about £>3 by an angle of f + A to generate an intermediate reference frame B defined by the unit 

III / 

vectors b l , b 2 , and b , (not shown), 2) the second rotation about b 2 by the dihedral angle T of the elastic axis that 
results in the intermediate reference frame c' defined by the unit vectors c\, c 2 , and (not shown), and 3) the third 
rotation about c, by an angle of K to result in the reference frame C. The transformation can be represented by a series 
of coordinate rotations as 


bi 

b 2 

b 3 


— sin A 

— cos A 

0 ' 


cosT 

0 

sinT 


’ 1 

0 

0 


Cl 

cos A 

— sin A 

0 


0 

1 

0 


0 

-1 

0 


Cl 

0 

0 

1 


-sinT 

0 

cosT 


0 

0 

-1 


. C 3 . 


— sin A cos r 

cos A 

sinAsinT 


Cl 

cos A cos r 

sin A 

— cos A sin r 


Cl 

— sinT 

0 

— cosT 


c 3 


The analysis can be repeated for the left wing. The reference frame D is aligned with the left wing’s elastic axis 
and is defined by the unit vectors d 1, d 2 , and d 3 . The B frame can be related to D through three successive rotations: 
1) the first rotation about —b 3 by an angle of f A to generate an intermediate reference frame B defined by the 

n n n rr 

unit vectors b l , b 2 , and b 2 (not shown), 2) the second rotation about b 2 by the dihedral angle T of the elastic axis that 

/ III 

results in the intermediate reference frame D defined by the unit vectors d j, d 2 , and d 2 (not shown), and 3) the third 

/ 

rotation about d x by an angle of n to result in the reference frame D. The relationship can be expressed as 
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B. Elastic Axis 

An analysis of the combined motion of the left wing is conducted in the present section, and the motion of the right 
wing is considered to be equivalent for symmetric flight. This analysis is equivalent to that in a previous study 4 5 and is 
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included for completeness. 

Let x represent the coordinate along the elastic axis of a wing running from root to tip. The wing pre-twist angle 
y(x) thus represents the incidence of the airfoil section at the corresponding elastic axis coordinate. A typical wing pre- 
twist varies from nose-up at the wing root to nose-down at the wing tip and is commonly referred to as a “wash-out” 
twist distribution. 

The internal structure of a wing is typically composed of a complex arrangement of load carrying spars and wing 
boxes that carry the stresses and strains introduced by aerodynamic forces and aeroelastic deflections. For this analysis, 
an equivalent beam approach is used which models the wing’s elastic behavior using equivalent stiffness properties. It 
is a common approach in analyzing aeroelastic deflections 7 and can be used to analyze high aspect ratio wings with 
good accuracy. The effect of wing curvature is ignored and straight beam theory is used to model the wing deflection. 
The axial or extensional deflection of a wing is also generally very small and is neglected. 



Figure 4. Left Wing Reference Frame 


Consider an airfoil section on the left wing as shown in Fig. 4 undergoing bending and torsional deflections. Let 
(x,y,z) be the coordinates of point Q on the wing airfoil section. The undeformed local airfoil coordinates of point Q 
are 
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where rj and c are the local airfoil coordinates, and y is the wing section pre-twist angle, positive nose-down. 12 The 
wing pre-twist is defined with respect to the elastic axis of the wing. 

Differentiating with respect to x gives 


y x 

r 

— siny 

— cosy 

V 


-zy 

Zx 

= 7 

cosy 

— siny 

_ 4 _ 


1 

1 


Let © be a torsional twist angle about the x-axis, positive nose-down, and let W and V be flapwise and chordwise 
bending deflections of point Q, respectively. Then, the rotation angle due to the elastic deformation can be expressed 
as 

<j>(x,t) = &di-W x d 2 + V x d 3 (5) 

where the subscript x denotes the partial derivatives of ©, W, and V with respect to x. 

Let (xi,yi,zi) be the coordinates of point Q on the airfoil in the reference frame D with aeroelastic deformation. 
The coordinates (xi,yi,zi) are computed using the small angle approximation as 
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Differentiating x\,y\, and z.\ with respect to x yields 
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Neglecting the transverse shear effect, the longitudinal strain is computed as 13 
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Ignoring the second-order terms and using the Taylor series expansion, sj * is approximated as 
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For a small wing twist angle y, (y ) 2 « 0. Then 
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The moments acting on the wing are then obtained as 13 
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where E is the Young’s modulus; G is the shear modulus; y is the derivative of the wing pre-twist angle; I yy , I yz , and 
I zz are the section area moments of inertia about the flapwise axis; J is the torsional constant; and B \ , /T, and B 3 are 
the bending-torsion coupling constants which are defined as 
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(15) 


The strain analysis shows that, for a pre-twisted wing, the bending deflections are coupled to the torsional deflection 
via the slope of the wing pre-twist angle. This coupling can be significant if the wash-out slope y is dominant as in 
highly twisted wings such as turbomachinery blades. 
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C. Coupled Bending-Torsion Equations 


Without considering chordwise bending of the wing, the equilibrium conditions for bending and torsion are expressed 
13 


as 
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where m x is the pitching moment per unit span about the elastic axis, /, is the lift force per unit span, and m y is the 
bending moment per unit span about the flapwise axis of the wing. 


1. Aerodynamic Forces and Moments 

Because the structural modeling is intended for use in a static aeroelasticity model, a steady-state aerodynamics model 
is used. Aerodynamic information can be obtained through vortex-lattice modeling to develop the forces and moments 
for coupled bending-torsion of a flexible wing. 



Figure 5. Airfoil Forces and Moments 


Neglecting the effect of downwash that is caused due to lift generation over a three-dimensional finite-wing, the 
sectional lift coefficient for an airfoil cross section, assuming linear aerodynamics, is as follows: 


c L (x) = c La (x)a c (x) 


( 18 ) 


where a c is known as the aeroelastic angle of attack and comprises of the rigid-body angle of attack and the contri- 
bution due to aeroelastic deformation. Note that a c is defined relative to the elastic axis of the wing and q, is the 
sectional lift coefficient in the wing reference frame. 

Let a be the aircraft’s rigid-body angle of attack and a e be the effect on the local angle of attack due to aeroelastic 
deformation at the aerodynamic center of the airfoil section. Note that these values a, a e are defined relative to the 
aircraft pitch axis, not the elastic axis. 
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where cl 0 is the zero angle of attack lift coefficient for the airfoil section. 

It is also important to note that the elastic contribution to the local aeroelastic angle of attack, a e , can be used to 
characterize an aeroelastic deformation. Given a deformation characterized by elastic axis twist © and vertical bending 
slope W x , the elastic contribution to the aeroelastic angle of attack can be calculated as 


a e {x ) = — ©(x)cosAcosr — W t (x)sinA 

where a e is about the aircraft pitch axis. 

The steady-state drag coefficient can be modeled by a parabolic drag polar as 
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where co 0 is the section parasite drag coefficient and k is the section drag polar parameter. This can be expressed in 
terms of rigid-body and aeroelastic contributions: 


c D {x) = c Dr (x) + c De (x) 
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The pitching moment coefficient about the aircraft pitch axis can be computed as 
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where e is the location of the aerodynamic center relative to the elastic axis defined in the streamwise direction 
perpendicular to the pitch axis, positive when the aerodynamic center is forward of the elastic axis, and c,„ ac is defined 
about the pitch axis, positive nose-up. 

The lift force, drag force, and pitching moment about the aircraft pitch axis are expressed as 
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where cos A takes into account the correction due to the elastic axis sweep, but is not needed in 
calculation since Eq. 27 is already about the pitch axis. 

The forces and moments in the local coordinate reference frame are obtained as 
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For a model with only flapwise bending and torsion considered, the beam deflection analysis is affected only by 
the terms /“, in " , and in" . The aerodynamic force and moment terms are thus considered to be 
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2. Propulsive Forces and Moments 

For wing-mounted engines, both the engine mass and thrust must be accounted for in the wing structural analysis. The 
propulsive force and moment vector are expressed in the reference frame D as 
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m e = r e x f e = (x e d\ - y e d 2 -z e d 2 ) xf e = 8(x-x e ) 
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where T is the engine thrust, m e is the engine mass, (x e , y e . z e ) is the coordinate of the engine thrust center such that y t , 
is positive forward of the elastic axis and z e is positive below the elastic axis, and 8 (x — x e ) is the Dirac delta function 
such that 
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3. Summary of Coupled Bending-Torsion Partial Differential Equations 

The force and moment equations due to aerodynamics and propulsive sources and inertial effects are formulated. 
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Inserting Eq. 13 and the force and moment terms into the governing equilibrium equations, Eqs. 16 and 17, the 
following equations which describe the coupled bending and torsion motion of the wing can be formulated: 
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IV. Finite-Element Modeling 


The development of the coupled bending-torsion partial differential equations allows for wing bending and tor- 
sional deflections to be solved. FEM, 8 a numerical technique that uses locally-defined basis functions to numerically 
approximate the solution of the governing partial differential equations, is used. The wing structure is discretized into 
n equally spaced one-dimensional beam elements, and FEM is applied. The bending and torsional deflections can be 
approximated as 

®(*,0 = E©i(*,0 (57) 

i= 1 

w(x ,t) = (58) 

i=i 

where i refers to the z'-th element. 

For each element, the bending and torsional deflections are approximated as 
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where the subscripts 1 and 2 denote values at nodes 1 and 2, and y/,(x), j = 1,2 and <j>k(x), k = 1,2, 3, 4 are the linear 
and Hermite polynomial shape functions 
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where x £ [0, 1] is the local coordinate and / = Ms the element length. 

The weak-form integral expressions of the coupled bending-torsion partial differential equations are obtained by 
multiplying the equations by Ng (x) and (x) and then integrating over the wing span. This yields 

V [ nJ^-{\gJ + EB\ (y') 2 l N 0 d t - EB 2 y N" v Wj\ dx = 
i = Qj 0 ax KL J ) 

n rl n rl 

^ / Ng (c m c) qooCos 2 Acdx+ ^ / Ng(—mge cg + mr^Ng 0,- — me cg N w Wj)dx 

Z4-/0 ,=1 7o 


+ L ^0 ['"<■ Cy« + )^0 0 < “ r=JC + 


z— 1 


At r 

^ Ng Ty e sin Ar +Tz e cos A — m e gy e + (Tx e cos A + Ty e sin A + m e gy e T)N w Wi 


i= 1 


(67) 
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£ ^ Nl d dxl (-EB 2 yN e 0 ; + EIyXwi)dx = 

« W / dc \ n 

Y Nj, ( cl -j - tan Ac I cos 2 A cdx + Y Nj(~PgA — pAN w Wi + pAe cg N g 9,)dx 

i=\J0 V dx ) , = | 

— m e N w Wi + m e y e NeOi + T sinAr — m e g+T cos A(N g 0i + j) + {T sinA + m e gT)N w Wi 




(=1 


rl _ 

+ £ N l ( Tx e cos A + Ty e sin A + m e gy e T) (N e 0, + y ) 


i=l 


(68) 


The expressions on the left hand sides can be integrated by parts upon enforcing the boundary conditions resulting 
j N lj~ x { [g7 + £Bi(7) 2 ] N g 9, - EB 2 y N'J Wi } dx = - j N g T { [GJ + EB^y) 2 ] N e 9 i -EB 2 yN'',w i }dx (69) 

rl rl 

Jo N l^- EB ^ N '^i+EIyyN' w Wi)dX= N7( 


, (—EB 2 y N g 9, + EI yy N w w,)dx = J N w l {—EB 2 y N e 0, + EI yy N w Wi)dx 

The elemental mass matrix, stiffness matrix, and force vector are then established as 


M; = / m 

r 2 k N J e Ne 

e cgNj)AI w 

dx + m e 

' (y 2 e+zl)N r e N 6 

y e N^N w 

Jo 

£cgN W N q 

nJn w 


- yeNjNe 

N? V N W 


(70) 


(71) 


Ki = 


1" 


GJ + EB\(y ) 2 

frp / t f'T Jt ~ 

N e T N e —EB 2 y NqN w 

Jo 


/ r/’T 

-EB 2 yNj 

N' e £/XX . 


dx 


0 


-T cos A NjN g - (Tx e cosA + Ty e smA + m e gy e r)Nj,N g 


f' 

/ m 8 

e c g Nl 

dx+ [ ( 

CjiiCNq 

+ 

0 

/ r T' 

1 

Jo 

. N w _ 

Jo y 

af 


c m c tan ANJ 


( Tx e cos A + Ty e sin A + m e gy e r)N^N w 
~(T sin A + m e gr)Nj,N' w 


cjoo cos 2 Acdx 


J X=X e 

(72) 


CmCtanAA^ 


q 0 c cos" A cdx 


+ 


-( Ty e smAT+Tz e cosA-m e gy e )Nl 
T sinAr — m e g+ T cos Ay+ (Tx e cos A + Ty e sinA + m e gy e T)y 


N 1 

1 ~ \ A > 


(73) 


(74) 


The globally assembled system is described by the matrix equation 

Mx e + Kx e = F 

where x e = 0! wi 9 2 w 2 w 2 0, w, w ; ... 6 n + 1 w n +i w n+1 

Equation 74 represents the governing equation for solving the structural deflection of a flexible wing given aero- 
dynamic and propulsive force and moment inputs. By setting x e = 0 , the equilibrium solution can be obtained through 
inverting the stiffness matrix and pre-multiplying the force matrix. This represents the static structural deflection of 
the wing based on the prescribed load input at that flight condition. 


x e = K~ l F 


(75) 


Information extracted from the wing solution state vector x e is used to characterize the wing’s structural deflection 
along the elastic axis of the wing © and W. It can also be used to calculate the elastic contribution to the aeroelastic 
angle of attack a e in Eq. 22. 
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V. Vortex-Lattice Aerodynamic Modeling 


Vorview is a computational tool used for aerodynamic modeling of aircraft configurations using vortex-lattice 
method. 14 Based on lifting line aerodynamic theory, Vorview provides a rapid method for estimating aerodynamic 
force and moment coefficients. Input vehicle geometries are discretized within Vorview into a series of panels which 
are then represented by placement of spanwise and chordwise locations of bound or horseshoe vortices. Vorview 
computes the vehicle aerodynamics in both the longitudinal and lateral directions independently, and these can be 
combined to produce the overall aerodynamic characteristics of the vehicle at any arbitrary angle of attack and angle 
of sideslip. 

Vorview is considered a medium fidelity tool, and limitations associated with vortex-lattice modeling in general 
apply to Vorview aerodynamic analysis. The drag prediction by Vorview is most reliable only for induced drag 
prediction due to the inviscid nature of any vortex-lattice method. Prediction of viscous drag due to boundary layer 
separation and wave drag due to shock-induced boundary layer separation are generally not conducted by vortex- 
lattice, and viscous drag must be estimated using other methods. 

In addition to force and moment analysis, Vorview can provide a rapid estimation of aerodynamic derivatives 
including dynamic derivatives due to angular rates. These aerodynamic stability and control derivatives are useful in 
analyzing the stability and handling characteristics of an aircraft configuration. Owing to the computationally efficient 
vortex-lattice method, aerodynamic derivatives can be estimated in Vorview fairly quickly. A flight dynamic model 
for a given vehicle configuration can be easily developed with Vorview, using the results from these stability and 
handling analyses. Vorview has been validated by both wind tunnel data 10 as well as the NASA Cart3D tool, 15 which 
is a high-fidelity inviscid (Euler) CFD analysis code targeted at analyzing aircraft performance in conceptual and 
preliminary aerodynamic design. In general, both Vorview and Cart3D seem to have similar predictive capabilities 
when compressibility is not a factor. 

Figure 6 illustrates an aerodynamic model of the GTM in Vorview. 



Figure 6. GTM Aircraft Model in Vorview 


In this study, Vorview will be used as the primary tool for conducting aerodynamic modeling for the aircraft 
configurations. Total aircraft characteristics as well as sectional data along the aircraft wing surfaces can be post- 
processed from Vorview. 


VI. Parasitic Drag Modeling 

Due to the inviscid nature of vortex-lattice drag estimates, a conceptual parasitic drag model will be utilized for 
development of the longitudinal trim model. This is necessary for realistic trim thrust values to be calculated, which 
must take into account viscous drag over the ESAC model’s airframe. 
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The skin friction model for all the surfaces will be done by calculating the arc lengths of the surfaces along the 
direction of flow and treating those lengths as flat plate surfaces. For airfoil surfaces, the arc lengths can be divided into 
top and bottom arc surfaces. For revolved surfaces, the model will calculate the skin friction coefficients by revolving 
the mean arc lengths. 

The Reynold’s number for an arc length c is calculated as 


Re c 


pVc 


(76) 


A critical Reynold’s number Re Xc 
calculate the critical length as 


is assumed to be Re Xc = 600,000. This critical Reynold’s number is used to 


= 


Re Xc p 

pV 


(77) 


The kinematic viscosity model is calculated in ^ by the Sutherland’s viscosity model: 


At = 0.3170 x 10~ 10 r2 


/ 734.7 \ 
\T + 216 ) 


(78) 


The revised skin friction coefficient for a single arc length is calculated as 


c f = 


1.328 x c 
c 

C f , laminar 


+ 0.072 Rep 2 - 0.072 — Rep~ 

c c c 

' v " 

C f , turbulent 


(79) 


The total skin friction coefficient for the aircraft components can then be formulated by integrating the skin friction 
coefficients either over the entire span of a wing, or the entire revolution of a surface. 

If pressure drag is to be conceptually estimated as well, then a form factor for the wing needs to be estimated. A 
form factor FF can can be defined for a wing surface as 16 

FF= + (“) +100 (") 4 ) (1-34M 018 (cosA,„) 0 - 28 ) (80) 

\ \X/C)m 'C/ \C / J 

where l - is the maximum thickness of the airfoil, (f.) ln is the location of the airfoil maximum thickness point, and A,„ 
is the sweep of the maximum-thickness line. 

For bodies of revolution, the form factor is evaluated as 16 


FF = 


60 / 
P + 400 


(81) 


where / is a fineness ratio given by / = , 1 and A max is the largest cross-sectional area perpendicular to the flow. 

V n^max 

The total parasitic drag coefficient including skin-friction and pressure drag due to separation and the supervelocity 
effect is calculated as 16 

C 0 o =^ Cf{FF ^ Q)Swet (82) 

S ref 

where all the components’ contributions to parasitic drag are summed up. The value S we , represents the wetted surface 
area of the component, and Q represents an interference factor. 

The parasite drag modeling method introduced here is also used to develop the wing’s spanwise parasitic drag 
coefficients used in Eq. 23 of the aeroelastic model. 


VII. Automated Geometry Modeling Tool 

An automated geometry generation tool is developed in Matlab and is used to close the loop between the structural 
and aerodynamic modeling needed to generate an aeroelastic model. The geometry generation tool uses structural 
deflection data that is computed by the FEM model and applies it to the undeformed aircraft wing geometry to reflect 
static aeroelastic deflections. The vehicle geometry modeler directly outputs a geometry input file that can be read by 
Vorview when computing an aeroelastic solution. 
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Figure 7. GTM Coordinate Systems 


Consider the reference frames in Fig. 7. The coordinate reference frame (xa^a^Za) defines the Body Station (BS), 
the Body Butt Line (BBL), and the Body Water Line (BWL) of the aircraft, respectively. The coordinate reference 
frame (xy . yy ,z.y) is the translated coordinate system attached to the nose of the aircraft such that xy = xb~ 13.25 ft, 
yv = >7;, and zv = Zb — 15.8333 ft. This reference frame is used by the vortex-lattice aerodynamic modeling tool. 
The aircraft body reference frame {xb^biZb) is the same B coordinate system defined earlier in Fig. 3 by the unit 
vectors b \ , b 2 , and b 2 . The B coordinate frame is attached to the aircraft center of gravity (CG) such that x K = xy—xy, 
y B =y v -yy, and zb = Zv —Zy, where (xy : yyzy) is the coordinate of the CG in the (xy,yyzy) reference frame. 11 

The vehicle geometry modeler has access to the outer mold line of the aircraft geometry. It is capable of applying 
geometric transformations onto the outer mold coordinates of the wing’s jig-shape to simulate aeroelastic deflection. 
Neglecting chordwise bending deflection and utilizing the coordinate system of the left wing developed earlier (coor- 
dinate frame D), the aeroelastic deflections in bending and torsion are expressed in a vector form as 

<$> = ®d x -W x d 2 (83) 


A r = — W sin W x d \ + W cos W x d 2 (84) 

The coordinate reference frame (x,y,z) of the left wing is related to the coordinate reference frame (xy . yy ,zy) by 
the relationship 


d 1 ' 


— sin A cos r 

— cos A sin r 

-sinT 




d 2 

= 

— cos A 

sin A 

0 


b 2 


d 2 


sin A sin y 

cos A sin r 

— cosT 


b 2 




— sin A cos r 

— cos A cos r 

— sinT 


-Vl 



= 

— cos A 

sinA 

0 


V 2 



sin A sin r 

cos A sin r 

— cosT 


-V3 


where (vi, v 2 , V3) are the unit vectors for the Vorview coordinate reference frame Gy , yy , z.y ') . 

Thus, the aeroelastic deflections result in a wing twist expressed as an incremental angle of attack Aa (positive 
nose-up), a horizontal deflection A yy (positive deflection towards wing tip), and a vertical deflection A zv (positive 
displacement upward): 

Aa = — ©cosAcosT — VF T sinA (86) 

A yy = —W sin cos AcosT — W cosVT r cosAsinr (87) 

Az v = — W sin sin r + W cos W x cos T (88) 

A coordinate transformation to account for wing aeroelastic deflections is performed by rotating a wing section 
about its elastic axis by the incremental angle of attack Aa and then translating the resultant coordinates by the 
horizontal deflection A)y and the vertical deflection Ay/. 

Note that the value of Aa is equivalent to the value of a e , or the local change in the angle of attack for a wing 
section due to aeroelastic deformation represented by Eq. 22. 
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VIII. Static Aeroelastic Model 


In a standard static aeroelastic model, it is understood that the modeling effort must take into account that aeroe- 
lastic structural deformations will affect aircraft aerodynamics, while changing the aerodynamics will thus change the 
structural deformations. In development of an aeroelastic model, it is crucial to include coupling between a struc- 
tural and aerodynamic model. Previous studies have analytically constructed fully coupled aeroelastic finite-element 
models that utilize rigid wing lift-curve slopes as an aerodynamic model. 4 - 5 

For this study, a static aeroelastic model is developed by closing the loop between a structural FEM model and a 
vortex-lattice aerodynamic model. For a model involving only flapwise bending and torsion, the aeroelastic deflection 
can be summarized by the quantities of aeroelastic elastic axis twist 0(x), aeroelastic vertical (flapwise) bending W(x), 
and aeroelastic vertical bending slope W x (x). These quantities are emphasized to be aeroelastic deflections, while the 
terms 0(x), W (x), W x (x) are considered structural deflection terms, which may or may not be the aeroelastic solution 
for a given flight condition. By closing a static aeroelastic loop as shown in Fig. 8, the structural deflections 0(x), 
W(x), W x (x) are expected to converge to 0(x), W(x), W x (x) as iterations are conducted. The structural and aeroelastic 
deformations can also be represented by the elastic contribution to the aeroelastic angle of attack in Eq. 22, or a e (x) 
and a e (x). 

The static aeroelastic code maps an input flight condition corresponding to an angle of attack a, Mach number 
M, engine thrust value T, and altitude h into the respective static aeroelastic deflection solution a flexible wing would 
experience at that flight condition. The following procedure is followed: 



Figure 8. Static Aeroelastic Model Concept 


1. Vortex -lattice modeling is conducted on an input geometry at the specified flight condition to determine the 
aircraft total aerodynamic quantities, as well as sectional coefficients cl(x), c mac (x), k(x), cl u (x) , and x ac (x) or 
the location of the section aerodynamic centers. 

2. The structural FEM model uses the sectional aerodynamic load inputs to calculate the wing’s structural deflec- 
tion 0(x) and W (x). 

3. The geometry generation tool converts 0(x) and W (x) into the series of deformations a e (x), Ayy(x), and A Zv{x), 
and generates a new aircraft geometry with the deformed wing. 

4. Steps 1-3 are repeated until a convergence criteria is met. 

This static aeroelastic model outputs 0(x), W’(x), W x (x), and tt e {x), and vortex-lattice can be used to determine the 
total aircraft coefficients such as Cl, Cp, and C m for the model with aeroelastic deflection. The aeroelastic convergence 
criteria is placed such that the aircraft Cl is evaluated at each iteration and the loop is exited if AC/J between iterations 
is below a certain tolerance. 
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IX. Longitudinal Trim Model 


The condition known as “trim” refers to the condition where an aircraft is in a state of constant velocity equilibrium. 
For longitudinal trim, this formulates a three degree-of-freedom system where the input variables are the aircraft’s 
angle of attack a , the elevator deflection 8 e (positive deflection is downward deflection), and the engine thrust T. 

The forces and moments are developed in the aircraft stability axes shown in Fig. 9. Only symmetric flight within 
a vertical plane of a non-rotating flat Earth is considered. 



w 



Figure 9. Aircraft Longitudinal Forces and Axes 

The body-fixed coordinate system is consistent with the previous usages of the B coordinate frame, while the 
stability or wind axes are represented by (xs,ys,zs) where x$ is in the direction of V a , the aircraft velocity. The axis 
yi, represents a vector parallel to the flat Earth surface and perpendicular to the aircraft vertical altitude. The value y g 
represents the aircraft glide path angle, and the pitch angle 9 p is defined as the angle between the aircraft forward body 
axis and the horizontal y/, direction. The aircraft is modeled to experience a lift force L, drag force D, weight W, and 
pitching moment M acting at the aircraft CG. The equilibrium equations in the stability axes are: 


^F Xs = 2T cos(a + e) — D — VTsin(y) = 0 

(89) 

^F zs = 2 T sin(a + e) +L — Wcos(y) = 0 

(90) 

My s = M + 2Tz e cose + 27’x e sine = 0 

(91) 


where the value T represents a single engine’s thrust value and hence a factor of 2 appears in the equations of motion. 
The value e is the engine mount angle, positive upwards. 

The equations are expanded using non-dimensional coefficients: 

fx s = 2T cos(a + e) - (C Dr + C De + C DSe 8 e )q«,Sref - W sin y g = 0 (92) 

fz s = 2 T sin (a + e) + (C Lr + C Le + C Lge 8 e )q 0 oS ref - W cos % = 0 (93) 

m ys = (C„, r +C me +C mSe 8 e )qooSrefC + 2Tz e cosE + 27’x f .sine = 0 (94) 

Let the static aeroelastic deformation be represented by a e = a e (ys , Ct, T), a function of the wing spanwise station 
ys, the aircraft angle of attack a, and the engine thrust T . For a model with wing flexibility, it must be taken into 
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account that the wing’s aeroelastic deformation will change the direction of the thrust vector. Let the coordinates 
(x es ■ y es ■ Z e . s ) represent the location of the engine thrust center in the aircraft stability axes such that y e$ represents the 
location of the engine thrust center along the spanwise axis, x es represents the distance of the engine thrust center in 
front of the aircraft center of gravity, and z e$ represents the distance of the engine thrust center below the aircraft center 
of gravity. The equations can then be expanded as 


fx , 5 — 2 T cos (a + e + <x e (y es )) — I Cd 0 + K (Clq +Cl u <x ) 


Cd v 


+2 K (C Lq + C La a) J C La a e dy s + k(^Jc Lci a e c/y s ^ + c o 5e 8 e j qooS ref - W sin y g (95) 


c De 


fzs — 2Tsin(a + e + oc e (y es )) + +Q, a o: + J CL a oc e dys+CL Se S ( ^JqooS re f — Wcos'Y g 


CLr 


CL e 


(96) 


v iy s — ^ (-'mu T c,„„ O' f- j (-'m a cc,.dys f 8^ q^S r( >fC 2T Ze cos ( € -j- (x r (y es ) ) T- 27’x e sin ( ft 4- cx { , (ye$ ) ) (97) 


Newton’s method will be used in an iterative fashion to solve the nonlinear system of equations. Newton’s method 
is well-established as one of the standard ways in which nonlinear systems of equations are solved and is dependent 
on calculation or approximation of the Jacobian for the system of equations. When the Jacobian is well-formulated, 
Newton’s method has a quadratic rate of convergence. 17 The basic equation for Newton’s method for a system of 
equations where the solution vector is x t is given by 


xf +1 =xM^r7(4) (98) 

When applied to the trim system, the solution state is given by x t = { a 8 e T } T - The Jacobian of the system 
of equations is represented as J, f is a vector / = { f x f :s m ys } T representing the values of the nonlinear system 
of equations, xf represents the k-th iterative solution, /(xf ) represents the values of the nonlinear system of equations 
at the k-th state, and [7(xf )] represents the Jacobian matrix evaluated at the A-th state. Through iteration, the values 
of xf +1 subsequently approach the solution such that / = { 0 0 0 } and the solution vector approaches the trim 

solution x, defined by a, 8 e , and T . 

Taking the partial derivatives of the aircraft longitudinal equations of motion allows us to populate the terms in a 
Jacobian matrix: 


7n — ^ —-2T sin(a + e + a t .(y e5 )) 

- (2 KCl a a + 2 KC La J C La a e dy s + 2K(C Lo +C La a + l) d f CL «^ dys j qooSref (99) 

h 2 =J§ L = - CD 8 e q°° S ref ( 10 °) 

-/13 = =2 cos (a + e + a e (y es )) — 2 T sin(a + e + a e (y es )) — ~ s ^ 

- (2K(C k) +C La cc+l) d $ CL ““ edys ^ qooSref (101) 

J 21 = ^ = 2Tcos(a + £ + a e (y es )) (l + + ( C La + ^ ref (102) 
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J ™=ji =c ^°° s -f 


d fzs 


J 23 — — 2sin(o: + £ + a e (y es )) + 2 cos (a + £ + cc e (y es )) 


da e (y es ) d f C La a e dy s 


dT 


dT 


qooS re f 


(103) 

(104) 


■A31 


dllly s 

da 


C ma qooSrefC+ d C "'a^edy s qooSrefC 

- 2 Tze S in ( £ + a e (y es ) ) ^ ^ ^ ^ + 2 Tx e cos ( e + a g (y e5 ) ) ^ ^ 

dfUy s 

■A' 2 = ^ ^ = C m g^C[ooSrefC 


(105) 

(106) 


■A33 = =2 (cos(e + a e (y es ))z es + sin (e + d ^ C '" a ^ ye ^ d} qoo S r efc 

+ 2 T (cos(e + a c (>' e5 ))x f > 5 - sin(e + a e (y es ))z es ) 

( d- 7 dx \ 

+ 2 T ( cos(£ + a e (y es ))-^|f +sin(e + a e (;y e5 ))^f J (107) 

The elastic terms in the Jacobian, which require running a static aeroelastic sensitivity study to a and T at each 
iteration, can be costly to calculate. Therefore, for the scope of this study, these terms are not used in approximation of 
the Jacobian. The value of 2KCj^a is also approximated as the value Co a , which is rapidly estimated by vortex-lattice 
code. As a trade-off in using these simplifications, the convergence rate of the trim algorithm is expected to decrease. 
The Jacobian for the system becomes: 


J = 


—2Tsin(a + s) — CD a q°aS re f — CD Se q°oS re f 2cos(a + e) 

2T cos(a + e) +CL a qooS re f Cl s q°oS re f 2sin(a + e) 

Cm a Q°°SrefC Cm^qooSrefC 2(cOs(£)7 (?( . 3in(£)x r ,<.) 


(108) 


The static aeroelastic framework in Fig. 8 is augmented with an additional step in which the Newton’s method 
trim procedure is added within the structural-aerodynamic loops. The resulting framework is represented in Fig. 10. 



Figure 10. Static Aeroelastic Longitudinal Trim Model 

Additional steps and modifications are added to the static aeroelastic approach to generate the aeroelastic trim 
procedure: 
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1 . Vortex -lattice modeling is conducted on an input geometry at an estimated trim flight condition to determine the 
aircraft total aerodynamic quantities and control derivatives ( Cl s ,C mg ,Co Se ) as well as sectional coefficients 
cl{x), c mac (x), k(x), c La (x), and x ac (x). 

2. The structural FEM model uses the sectional aerodynamic load inputs to calculate the wing’s structural deflec- 
tion 0(x) and W (x). 

3. The geometry generation tool converts 0(x) and W (x) into the series of deformations a e (x), A yy(x), and A Zv{x), 
and generates a new aircraft geometry with the deformed wing. 

4. Newton’s method is used to determine a new trim flight condition and inputs a, S e , and T . 

5. Steps 1-4 are repeated until a convergence criteria is met. 

The outputs of the static aeroelastic trim code are the flexible aircraft trim inputs a, 8 e and f as well as the aeroelastic 
deflections 0(x), W(x), W x (x), and a e (x). 

The convergence of the trim code is evaluated based upon the non-dimensionalized values of the forces and moment 
in the stability axes. Let J3i, /3?, and J83 be user defined parameters. The following three equilibrium criteria must be 
satisfied for a trim solution to be determined: 


qooO re f 

(109) 

-4- ^ 1&1 

Qoa^ref 

(HO) 

qooS re fC 

(111) 


A convergence criteria is also placed on the aircraft quantity A('C/,„ + Ci f: ) such that its value between iterations 
must also be below a certain threshold. This differs from the static aeroelastic model convergence criteria on |ACl| 
because the aircraft angle of attack a is not held constant in trim model. This change in convergence criteria helps to 
ensure that aeroelastic shape is converged as the trim solution is determined. 

The framework in Fig. 10 represents a simplistic aeroelastic trim procedure where Newton’s method is added di- 
rectly into the structural-aerodynamic iterations. Other variations of a trim algorithm could use structural-aerodynamic 
coupling as an inner loop while Newton’s method is wrapped around the static aeroelastic model, thus converging the 
aeroelastic shape of the wing at each flight condition prior to conducting a Newton’s method iteration. These ap- 
proaches are not investigated in this work. 

X. Simulations and Results 

The following sections summarize the results from simulating the static aeroelastic and aeroelastic longitudinal 
trim models for different aircraft models. Flexible wing structural and aeroelastic deflections are presented, and the 
previous used conventions are maintained: 

0(x) represents elastic axis twist, positive nose-down. 

W(x) represents vertical or flapwise bending, positive upwards. 

W x {x) represents vertical or flapwise bending slope, positive corresponding with positive upwards deflection. 

a e (x) represents elastic contribution to the angle of attack or twist about the pitch axis, positive nose-up. 

Values represented with a bar such as © and W represent aeroelastic deflections computed through convergence of 
structural-aerodynamic iterations. 

A. Idealized Wing Alone Model, Static Aeroelastic Analysis 

The static aeroelastic code is implemented on the Idealized Wing Alone model as a validation study. The Idealized 
Wing Alone model is simulated at the 80% fuel cruise case of the GTM, corresponding to M = 0.8, h = 30,000 ft, 
and W = 190,000 lbs. Representative mass and stiffness values based on the GTM data are used to generate the 
beam model for the Idealized Wing Alone, and the empty weight of the wing is simulated to be W w fi — 17,250 lbs. 
The weight of the engines are simulated on the structural beam model as point loads, but the thrust effect on the 
structural model are ignored. The aerodynamic model of the Idealized Wing Alone does not include the engine as a 
simplification. 


19 of 29 


American Institute of Aeronautics and Astronautics 


Lift curves and drag polars for the Idealized Wing Alone model are generated for different cases and are shown 
in Fig. 11. The drag polars for the Idealized Wing Alone do not include parasitic drag and only show the vortex- 
lattice computed drag. In the cases where the the wing is considered “rigid”, the wing’s structural deflections do not 
have any coupling to the aerodynamic modeling. Thus, for the cases labeled “rigid”, the aerodynamics on the model 
are equivalent to the loads that a rigid planform would experience. The structural deflection results are simply the 
deflection that would be expected based on the rigid wing loads. These differ from the coupled aeroelastic solutions 
corresponding to the results labeled as “flexible”. 




Figure 11. Lift Curves and Drag Polars for Idealized Wing Alone Model 

Figure 12 shows the deflection results for the flight condition at M = 0.8, for the Idealized Wing Alone model. The 
rigid models’ tip deflection W n[) are presented along with the flexible models’ aeroelastic tip deflections W Iip . The tip 
elastic twist values © r , p and 0,, p are also presented for the rigid and flexible models. 
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Figure 12. Structural and Aeroelastic Deflections for Idealized Wing Alone Model 

Based on the aircraft weight of W = 190,000 lbs, a trim C) for a 1-g load is determined for the Idealized Wing 
Alone to be Cl = 0.346. Because the Idealized Wing Alone is not modeled with a horizontal tail or with engine 
thrust, a simple angle of attack a trim can be conducted based off the lift curves developed in Fig. 11. The static 
aeroelastic deflections (or the structural deflections for the rigid wing case) are determined and compared against 
NASTRAN static aeroelastic results provided by Boeing Research & Technology. 2 The results for the rigid wing 
structural deflection solutions are presented in Table 1 . In addition to the wing tip deflection W tlp and wing tip elastic 
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axis twist 0 /jp , the value of the elastic contribution to angle of attack 0 . t , ip is also provided representing the amount 
of twist of the wing tip about the pitch axis. 


Planform Only 

Wing Alone Model 

NASTRAN 

Absolute 

Difference 

a, deg 

3.20 

3.15 

0.05 

W tip , in 

24.34 

22.42 

1.92 

©tip , deg 

-1.43 

-1.12 

0.31 

a eJip , deg 

-0.25 

-0.36 

0.11 


With Twist 

Wing Alone Model 

NASTRAN 

Absolute 

Difference 

a, deg 

2.33 

2.27 

0.06 

Wtip, in 

16.97 

15.12 

1.85 

©tip, deg 

-1.06 

-0.80 

0.26 

OLeJipi deg 

-0.07 

-0.16 

0.09 


With Camber 

Wing Alone Model 

NASTRAN 

Absolute 

Difference 

a , deg 

1.28 

1.16 

0.21 

Wu p , in 

26.50 

28.28 

1.78 

©tip, deg 

-0.05 

0.19 

0.24 

oce.tip, deg 

-1.64 

-1.94 

0.30 


Combined 

Wing Alone Model 

NASTRAN 

Absolute 

Difference 

a, deg 

0.41 

0.29 

0.12 

W tip , in 

19.05 

20.99 

1.94 

@ti P , deg 

0.32 

0.51 

0.19 

a e ,ti P , deg 

-1.46 

-1.74 

0.28 


Table 1. Structural Deflection (Rigid Wing) Results for Idealized Wing Alone, h = 30,000 ft 

The structural deflection results for the Idealized Wing Alone model when only the rigid wing loads are considered 
are in good agreement between the aeroelastic model and the NASTRAN results. The percent differences in tip 
deflection between the Idealized Wing Alone model and the NASTRAN results are 8.56%, 12.2%, 6.29%, and 9.24% 
respectively for the planform only model, planform and twist, planform and camber, and combined models. Wing tip 
twist agrees where the absolute differences of & t i P < |0.31°| and a eJip < |0.30°| are observed for elastic axis twist and 
pitch axis twist, respectively. 

The effects of adding twist, camber, and twist and camber on the aeroelastic model with no coupling serve as 
confirmation of expected aerodynamic behavior. Because the angle of attack a is adjusted to maintain the total load 
of W = 190,000 lbs, the effect of twist, camber, and both combined acts only to redistribute the lift load along the 
planform of the wing. With just the planform, the angle of attack to maintain Cl = 0.346 is a = 3.20°. With twist 
enabled, the angle of attack decreases to a = 2.33°. The decrease in angle of attack is expected due to the positive 
nose-up twist of the jig-shape at the root which increases the net lift on the planform. The wash-out of the wing twist 
shifts the lift distribution towards the root, and thus, a decrease in the tip deflection W r , p and more nose-down tip twist 
© r;p are observed. The introduction of camber also significantly affects the angle of attack for Cl = 0.346. It is known 
from aerodynamics that camber significantly increases the lift on airfoil sections and also introduces moment about 
the airfoil aerodynamic center. Because the cambered wing is able to generate much more lift than the flat planform, 
the angle of attack needed to maintain the 1-g lift load is only a = 1.28°, much less than the planform only. The 
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tip twist ® t ip also becomes more nose-down. This is expected, as typical cambered aircraft wings have negative c mac 
which twists the airfoil sections nose-down. 

The structural deflection results without aeroelastic coupling are recomputed for an altitude of h = 35,000 ft {Cl = 
0.437) and the results are summarized in Table 2. 


Planform Only 

Wing Alone Model 

NASTRAN 

Absolute 

Difference 

a, deg 

4.04 

3.98 

0.06 

W tip , in 

24.34 

22.42 

1.92 

©r i P , deg 

-1.43 

-1.12 

0.31 

Otejipi deg 

-0.25 

-0.36 

0.11 


With Twist 

Wing Alone Model 

NASTRAN 

Absolute 

Difference 

a, deg 

3.17 

3.09 

0.08 

Wup, in 

18.53 

16.63 

1.90 

@ti P , deg 

-1.14 

-0.86 

0.28 

a e jip, deg 

-0.10 

-0.20 

0.10 


With Camber 

Wing Alone Model 

NASTRAN 

Absolute 

Difference 

a, deg 

2.12 

1.99 

0.13 

Wup, in 

26.03 

27.07 

1.04 

©tip, deg 

-0.34 

-0.08 

0.26 

deg 

-1.35 

-1.61 

0.26 


Combined 

Wing Alone Model 

NASTRAN 

Absolute 

Difference 

a, deg 

1.25 

1.10 

0.15 

W ti p, in 

20.18 

21.28 

1.10 

©tip, deg 

-0.04 

0.18 

0.22 

(%e,tipi deg 

-1.21 

-1.46 

0.25 


Table 2. Structural Deflection (Rigid Wing) Results for Idealized Wing Alone, h = 35,000 ft 

Good agreement still exists at the increased altitude of h = 35,000 ft between the developed aeroelastic model 
and NASTRAN results for the Idealized Wing Alone model when only the rigid wing loads are considered. Percent 
differences of 8.56%, 11.4%, 3.84%, and 5.26% respectively for the planform only model, planform and twist, plan- 
form and camber, and combined models are observed between the wing tip deflection Wu p results with the developed 
aeroelastic model and NASTRAN. Absolute differences of ©tip ^ |0.3 1° | and a e , t i P < |0.26°| are observed for wing 
tip elastic axis twist and pitch axis twist, respectively. 

The aerodynamic trends observed for the Idealized Wing Alone model with the rigid wing at h = 30, 000 ft are 
consistent with the results at h = 35,000 ft. The addition of twist, camber, and twist and camber to the planform 
increases the lift generated by the model and subsequently reduces the angle of attack a needed to maintain a load of 
C L = 0.346. 
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The results when coupling is enabled and the framework is true to Fig. 8 is presented in Table 3 for Cl = 0.346 at 
h = 30. 000 ft. 


Planform Only 

Wing Alone Model 

NASTRAN 

Absolute 

Difference 

a, deg 

3.44 

3.41 

0.03 

W tip , in 

23.13 

21.28 

1.85 

®ti P , deg 

-1.39 

-1.08 

0.31 

&e,tip’> deg 

-0.21 

-0.33 

0.12 


With Twist 

Wing Alone Model 

NASTRAN 

Absolute 

Difference 

a, deg 

2.47 

2.42 

0.05 

Wu p , in 

16.43 

14.62 

1.81 

&u P , deg 

-1.04 

-0.78 

0.26 

OC e ,t ip, deg 

-0.05 

-0.15 

0.10 


With Camber 

Wing Alone Model 

NASTRAN 

Absolute 

Difference 

a, deg 

1.84 

1.88 

0.04 

Wup, in 

22.07 

23.43 

1.36 

®ti P , deg 

0.12 

0.38 

0.26 

oc e ,ti P , deg 

-1.49 

-1.79 

0.30 


Combined 

Wing Alone Model 

NASTRAN 

Absolute 

Difference 

a, deg 

0.87 

0.90 

0.03 

W tip , in 

15.34 

16.78 

1.44 

deg 

0.49 

0.68 

0.19 

&e,tip-> deg 

-1.35 

-1.62 

0.27 


Table 3. Aeroelastic Deflection (Flexible Wing) Results for Idealized Wing Alone, h = 30,000 ft 

The good agreement between the aeroelastic Idealized Wing Alone model and the NASTRAN aeroelastic model 
results remains when structural-aerodynamic coupling is enabled. Percent differences of 8.69%, 12.3%, 5.80%, and 
8.58% between the wing tip deflection W t i P results of the developed aeroelastic model and the NASTRAN model are 
observed for the planform only, planform and twist, planform and camber, and combined models, respectively. For 
wing tip twist, absolute differences of ©„ p < |0.31°| and a eJ i P < |0.30°| are observed for wing tip elastic axis twist 
and pitch axis twist, respectively. 

The results in Table 3 can also be compared against the results in Table 1 to illustrate the effect of adding in 
aeroelastic coupling on the model. It is observed that the angle of attacks a for the flexible models at Cl = 0.346 
are higher than that of the rigid models. This is due to the fact that the aeroelastic deflection on the wings tends to 
cause flexible wings to twist downwards with a negative a e . This reduces the amount of lift generated by the wing 
at a particular angle of attack a. Thus, in order to maintain the same load as a rigid model, a flexible model needs 
to operate at a higher angle of attack a. For example, in the combined Idealized Wing Alone model where twist 
and camber are both considered, a higher angle of attack a = 0.86° is required for the flexible model while a lower 
angle of attack a = 0.41° is adequate for the rigid model for Cl = 0.346. This phenomenon is also illustrated in 
Fig. 11, where the flexible wing lift curves are shifted below that of the rigid wing lift curves. The exception is for 
the uncambered models at low angles of attacks, where the elastic contribution to the aeroelastic angle of attack a e 
can actually be positive, twisting the wing nose-up. This is due to the fact that the flapwise bending slope W x , which 
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generally contributes to the nose-down elastic twist of the wing, is very low at these flight conditions. 

The results of the static aeroelastic model for the increased altitude h = 35,000 ft or Cl = 0.437 flight condition 
are also compared against NASTRAN results. 


Planform Only 

Wing Alone Model 

NASTRAN 

Absolute 

Difference 

a, deg 

4.29 

4.23 

0.06 

W tip , in 

23.36 

21.50 

1.86 

®np, deg 

-1.40 

-1.09 

0.31 

OC e ,tip , deg 

-0.21 

-0.33 

0.21 


With Twist 

Wing Alone Model 

NASTRAN 

Absolute 

Difference 

a, deg 

3.34 

3.27 

0.07 

W ti p, in 

17.98 

16.13 

1.85 

©r;>, deg 

-1.12 

-0.85 

0.27 

u e ,tip, deg 

-0.09 

-0.19 

0.10 


With Camber 

Wing Alone Model 

NASTRAN 

Absolute 

Difference 

a, deg 

2.64 

2.62 

0.02 

Wnp, in 

22.90 

23.77 

0.87 

®i,p, deg 

-0.22 

0.04 

0.26 

Otejipi deg 

-1.24 

-1.51 

0.27 


Combined 

Wing Alone Model 

NASTRAN 

Absolute 

Difference 

a, deg 

1.67 

1.66 

0.01 

W n p, in 

17.57 

18.40 

0.83 

©r i P , deg 

0.07 

0.29 

0.22 

a e ,,ip, deg 

-1.13 

-1.37 

0.24 


Table 4. Aeroelastic Deflection (Flexible Wing) Results for Idealized Wing Alone, h = 35,000 ft 

Good agreement continues to exist between the results of the developed aeroelastic model and that of the NAS- 
TRAN generated results. The observed relative percent differences between aeroelastic wing tip deflections are 8.65%, 
11.5%, 3.66%, and 4.51% respectively for the planform only model, planform and twist, planform and camber. Ab- 
solute differences in aeroelastic wing tip twist about the elastic axis are observed to be |A0„- p | < 0.31°, and absolute 
differences in the elastic contribution to aeroelastic angle of attack are observed to be |Aa e r , p | < 0.26°. 

Based on the results in Tables 1-4, it is observed that the elastic axis twist of the developed aeroelastic model 
consistently experiences more nose-up twist than that of the NASTRAN results. Nonetheless, the good agreement for 
all the test cases produces confidence in the aeroelastic framework developed and its predictive capability for flexible 
wing aircraft. 
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B. ESAC Model, Static Aeroelastic and Longitudinal Trim Analysis 

With preliminary validation of the aeroelastic modeling completed, the static aeroelastic framework is extended to 
that of the ESAC. The ESAC model is representative of a full aircraft configuration with engine nacelles, pylons, tail 
surfaces, and fuselage present in the modeling. 

Two independent and separate cruise conditions are simulated with the ESAC, representing different possible de- 
sign cruise condition candidates. The weight model is also adjusted based on the different candidate cruise conditions. 

• A first cruise condition is simulated at M = 0.8, h = 30,000 ft, where the 80% fuel case of the aircraft is 
W = 190,000 lbs and the empty wing mass is W W) o = 17,250 lbs. The design lift coefficient in this case is 
C L = 0.346. 

• A second cruise condition is simulated at M = 0.8, h = 36,000 ft, where the 80% fuel case of the aircraft is 
W = 210,000 lbs and the empty wing mass is W w ,o = 13,000 lbs. The design lift coefficient for this case is 
C L = 0.510. 

For the cruise flight conditions, the parasitic drag coefficient is built up using a critical Reynold’s number of Re c = 
600,000. The results are presented in Table 5. 


h = 30,000 ft 

Component 

Cf 

C Do 

Wings 

0.0037 

0.0054 

Fuselage 

0.0042 

0.0047 

Horizontal Tail 

0.0010 

0.0015 

Vertical Tail 

0.0008 

0.0011 

Engine Nacelles+Pylons 

0.0004 

0.0010 

Total 

0.0102 

0.0137 


h = 36,000 ft 

Component 

Cf 

Cd 0 

Wings 

0.0038 

0.0056 

Fuselage 

0.0044 

0.0049 

Horizontal Tail 

0.0011 

0.0015 

Vertical Tail 

0.0008 

0.0012 

Engine Nacelles+Pylons 

0.0005 

0.0011 

Total 

0.0106 

0.0142 


Table 5. Parasite Drag Build-Up At Cruise Conditions 

Initially neglecting thrust and its effect on the aeroelastic deformation, the lift curve and drag polars for the ESAC 
are generated and shown in Fig. 13. Note that now the drag polar includes the conceptually estimated parasitic drag 
coefficient. 




Figure 13. Lift Curves and Drag Polars for ESAC Model, No Thrust or Control Deflection 

The tip deflection, W t i P for the rigid model and W lt/I for the flexible model, and tip elastic axis twist, © ri/7 for the 
rigid model and © f;p for the flexible model, at the two different flight conditions are shown in Fig. 14. 
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Figure 14. Structural and Aeroelastic Deflections for ESAC Model, No Thrust or Control Deflection 

The plots in Fig. 13 and 14 are consistent with expected behavior. The lift curve of the flexible wing models 
shifts below the rigid wing models due to the aeroelastic deformation’s effect on twisting the flexible wing nose down, 
effectively decreasing the local angle of attack along the wing span. The effect of coupling structural deflection with 
aerodynamic vortex-lattice also resulted in effective stiffening in bending and softening in torsion. This is observed by 
the decrease in flapwise deflection in comparing the flexible wing models to the rigid wing models and the increase in 
nose-down twist. 

The longitudinal trim framework is then applied to the ESAC model configuration, and control deflections of the 
elevator are enabled as well as thrust from the two engines. Deflection of the elevator surface is assumed not to affect 
the aerodynamics over the wing, but the thrust effect on the wing is considered. The trim algorithm framework can be 
readily applied to both the rigid and flexible wing models, and the rigid wing model is analyzed first for comparison. 
The convergence criteria for the trim algorithm is set such that /3i = /fe = 0.001 in Eqs. 109 and 1 10 and /3 3 = 0.00005 
inEq. 111 . 

The static trim results for the two cruise conditions when the ESAC’s wings are considered rigid are presented 
in Table 6 . Once again, the deflection results presented for the rigid wing solution represent structural deflections 
assuming that the deformation of the wings does not affect the aerodynamics. 


Cruise Condition #1 Design Cl = 0.346 

a, deg 

2.108 

8 e , deg 

-5.075 

T, lbf (single engine) 

5821 

L/D 

14.8 

W tip , in 

28.082 

®n P , deg 

0.459 

Ote,tip-> deg 

-2.209 


Cruise Condition #2 Design Cl = 0.510 

a, deg 

3.640 

S e , deg 

-6.254 

T, lbf (single engine) 

6249 

L/D 

14.8 

W,i P , in 

30.810 

©r ip, deg 

-0.125 

(%e,tip) deg 

1.841 


Table 6. Static Trim Results For Rigid Wing ESAC 

The rigid wing ESAC model’s trim solution is very easily determined by the trim code, and a total of three iterations 
of Newton’s method were conducted from a user specified estimated initial point of a = 2.0° for the first cruise 
condition of design Cl = 0.346 and the initial point of a = 3.0° for the second cruise condition. The pilot input trim 
states at each of these iterations are plotted and shown in Fig. 15. 
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| Cruise Condition #1 , Design C L = 0.346j 


| Cruise Condition #2, Design C L = 0.51 



4 



8 3.5 




Iteration Number 


Iteration Number 




Iteration Number 


Iteration Number 


Iteration Number 



Iteration Number 


Figure 15. Static Trim Iteration Data for Rigid Wing ESAC 

The longitudinal trim code is then run on the ESAC with fully flexible wings. The trim solution is determined, and 
the results are summarized in Table 7. 


Cruise Condition #1 

Design Cl = 0.346 

a, deg 

3.747 

S e , deg 

-6.798 

T, lbf (single engine) 

7054 

L/D 

11.4 

W t ip, in 

24.645 

®ti P , deg 

0.643 

(%e,tip-> deg 

-2.159 


Cruise Condition #2 

Design Cl = 0.510 

a, deg 

5.160 

S e , deg 

-7.880 

T, lbf (single engine) 

7285 

L/D 

11.9 

W tip , in 

28.469 

©tip, deg 

0.014 

a e ,ti P , deg 

-1.817 


Table 7. Static Trim Results For Flexible Wing ESAC 

Due to the flexibility of the wing, the trim code takes longer to converge, with 20 iterations for the design Cl = 
0.346, and 14 iterations for the design Cl = 0.510. The pilot input trim states at each of these iterations are plotted in 
Fig. 16. As expected, the trim angle of attacks a for the flexible models are higher than that of the rigid models due 
to the effect of the aeroelastic bending and torsion on the aeroelastic angle of attack. 

Because the flexible trim case also involves determining the coupled aeroelastic wing deformation, plots of the 
evolution of the tip deformation are shown in Fig. 17, where W t j P is plotted in feet and 0,, ;) is plotted in degrees. It is 
observed that the aeroelastic deformation of the model converges alongside the solution of the static trim inputs. 


27 of 29 


American Institute of Aeronautics and Astronautics 


Cruise Condition #1 , Design C L = 0.346 
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Cruise Condition #2, Design C L = 0.510 
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Figure 16. Static Trim Iteration Data for Flexible Wing ESAC 


Cruise Condition #1, Design C L = 0.346 


Cruise Condition #2, Design C L = 0.510 
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Figure 17. Static Trim Deflection Iteration Data for Flexible Wing ESAC 


XI. Conclusions 

This study presents the development of a static aeroelastic model and a static three degree-of-freedom longitudinal 
aeroelastic trim model capable of analyzing flexible wing aircraft. The aeroelastic model is developed by coupling 
equivalent beam wing structural models in the loop with vortex-lattice aerodynamic modeling. The static aeroelastic 
code is capable of generating flexible wing aircraft configurations and characterizing aeroelastic deformations. Utiliz- 
ing converged aeroelastic geometry, flexible aircraft lift curves and drag polars can be generated exploring the effect of 
flexibility enabled models. The coupled finite-element vortex-lattice aeroelastic model is validated against NASTRAN 
aeroelasticity results for a simplified Idealized Wing Alone model for various test cases, including systematically ana- 
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lyzing the model planform only, with only added twist, with only added camber, and both added. The static aeroelastic 
model is extended to develop an aeroelastic three degree-of-freedom longitudinal trim model for a full aircraft config- 
uration with horizontal tail, elevator, and engines. The trim algorithm uses Newton’s method which is inserted in the 
static aeroelastic approach to iterate pilot trim inputs and to solve the flight condition until convergence is achieved. 
This framework is readily applied for the flexible wing models investigated in this study, but could be improved in the 
future for highly flexible wings. Future study of an aeroelastic trim model where Newton’s method is wrapped around 
an inner loop static aeroelastic model can be investigated. In addition, different approaches for estimating the Jacobian 
and variations to Newton’s method can possibly increase the convergence rate of the longitudinal trim model. 
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